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Abstract 

in. 

£NJ , We study the nonexpansivity of reflection mappings in geodesic spaces and apply our find- 

ings to the averaged alternating reflection algorithm employed in solving the convex feasibility 
problem for two sets in a nonlinear context. We show that weak convergence results from Hilbert 
spaces find natural counterparts in spaces of constant curvature. Moreover, in this particular 
' setting, one obtains strong convergence thereof. 

1 Introduction 

The convex feasibility problem for two sets consists in finding a point in the intersection of two 
nonempty closed and convex sets provided such a point exists. This problem finds remarkable 
qq 1 applications in applied mathematics and various branches of engineering (see, for example, [3, 11, 

I/"") ■ 28, 10]) which motivated many researchers to focus on methods of solving this problem. In Hilbert 

spaces there exists a wide range of algorithms that use metric projections on the sets in order to 
obtain sequences of points that converge weakly or in norm (under more restrictive conditions) to a 
solution of this problem. One of the most famous algorithms is the alternating projection method 
£NJ . which was developed by von Neumann [30] and was recently adapted to the setting of CAT(O) 

spaces by Bacak, Searston and Sims in [2]. 

Another class of algorithms considered in this respect bases on reflections instead of projections. 
Given a nonempty closed and convex subset A of a Hilbert space H, the reflection of a point x G H 
with respect to A is the image of x by the reflection mapping Ra = 2Pa — I, where Pa stands 
for the metric projection onto A and / is the identity mapping. In this work we focus on the 
averaged alternating reflection (AAR) method employed in solving the convex feasibility problem 
for two sets. Suppose A and B are two nonempty closed and convex subsets of a Hilbert space with 
nonempty intersection. The AAR method generates the following sequence for a starting point 

xq E H: x n = T n xo, where T = — — ^ A ^ B , This algorithmic scheme was studied by Bauschke, 

Combettes and Luke in [4, 5] not only in connection with the convex feasibility problem, but also 
for finding a best approximation pair of the sets A and B in case their intersection is empty and 
such a pair exists. The AAR method was later modified in [6] in order to solve the problem of 
finding the projection of a point onto the intersection of two closed and convex sets. In fact, for the 
convex feasibility problem, this algorithm is a special case of one described by Lions and Mercier 
in [25]. One obtains weak convergence of the sequence (x n ) to a fixed point of the mapping T and 
the projection of this point onto the set B lies in the intersection of A and B. 

Here we are interested in studying the AAR method in geodesic spaces. However, several 
difficulties arise when considering this algorithm in a nonlinear setting. First of all one needs to 
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find an appropriate definition for the reflection mapping. In order to guarantee the existence of 
the reflection of a point in the space, we consider spaces with the geodesic extension property. 
Moreover, the reflection mapping is not always unique. A second difficulty consists in guaranteeing 
certain properties of this mapping. In Hilbert spaces, the proof of the convergence of the AAR 
method relies on the nonexpansivity of the reflection mapping which yields the firm nonexpansivity 
of the mapping T. 

In this paper we prove that the reflection mapping is nonexpansive in spaces of constant cur- 
vature and justify why it fails to be nonexpansive in the broad setting of CAT(O) spaces. We also 
analyze the behavior of reflection mappings in slightly more general settings, namely gluings of 
model spaces. Likewise, we study the convergence of the AAR method in spaces of constant curva- 
ture proving strong convergence thereof. Furthermore, we include a rate of asymptotic regularity 
for the AAR method. 

This work is partly motivated by a communication of Ian Searston given during the 10th 
International Conference on Fixed Point Theory and its Applications where the problem of studying 
the nonexpansivity of the reflection mapping was raised. 

2 Preliminaries 

A metric space (X, d) is said to be a (uniquely) geodesic space if every two points x and y of X 
are joined by a (unique) geodesic, i.e, a map c: [0,/] C R — > X such that c(0) = x, c(l) = y, and 
d(c(t),c(t')) = \t — t'\ for all t, t' £ [0,1]. The image c([0, Z]) of a geodesic forms a geodesic segment 
which joins x and y and is not necessarily unique. If no confusion arises, we use [x, y] to denote 
a geodesic segment joining x and y. A point z in X belongs to a geodesic segment [x, y] if and 
only if there exists t £ [0, 1] such that d(x,z) = td(x,y) and d(y,z) = (1 — t)d(x,y) and we write 
z = (1 — t)x + ty for simplicity. Notice that this point may not be unique. A subset of X is said to 
be convex if it contains any geodesic segment that joins every two points of it. A geodesic triangle 
A(x,y,z) consists of three points x,y,z £ X (the vertices of A) and three geodesic segments 
joining each pair of vertices (the edges of A). A geodesic line in A is a subset of X isometric to 
R. A geodesic space has the geodesic extension property if each geodesic segment is contained in a 
geodesic line. More on geodesic metric spaces can be found for instance in [7, 27]. 
The metric d: X x X — > R is said to be convex if for any x,y,z £ X one has 

d(x, (1 - t)y + tz) < (1 - t)d(x, y) + td(x, z) for all t £ [0, 1]. 

A geodesic space (X,d) is Busemann convex (introduced in [9]) if given any pair of geodesies 
c\ : [0, l\] — > X and c 2 : [0, 12] — > X one has 

d{c 1 {th),c 2 {tl 2 )) < (1 - t)d(c 1 (0),c 2 (0))+td(c 1 (l 1 ),c 2 (h)) for all t G [0,1]. 

It is well-known that Busemann convex spaces are uniquely geodesic and with convex metric. 

A very important class of geodesic metric spaces are CAT(fc) spaces (where k G R), that is, 
metric spaces of curvature uniformly bounded above by k in the sense of Gromov. CAT(fc) spaces 
are defined in terms of comparisons with the model spaces MJ}, which are the complete, simply 
connected, Riemannian ra-manifolds of constant sectional curvature k. Since these model spaces 
are of essential importance in this work we give their definition directly as metric spaces and recall 
some of their properties. For a thorough treatment of such spaces and related topics the reader 
can check [7, 17]. 

The n-dimensional sphere § n is the set {x £ R n+1 : (x \ x) = l}, where (• | •) is the Euclidean 
scalar product. Define d : S n x S n — > R by assigning to each (x,y) £ S n x S n the unique number 
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d(x,y) G [0, 7r] such that cosd(x,y) = (x \ y). Then (§ n ,d) is a metric space called the spherical 
space. This is a geodesic space and if d(x, y) < ir then there is a unique geodesic joining x and 
y. Also, balls of radius smaller than ir/2 are convex. The spherical law of cosines states that in a 
spherical triangle with vertices x,y,z G § n and 7 the spherical angle between the geodesic segments 
[x, y] and [x, z] we have 

cos d(y, z) = cos d(x, y) cos d(x, z) + sin d(x, y) sin d(x, z) cos 7. 

For u,v £ R n+1 , consider the quadratic form given by {u \ v) = —u n+ iv n+ i + ^27=1 u i y i- The 
hyperbolic n-space HP is the set {u = (ui,U2, ...,u n +i) G R n+1 : (u | u) = — l,u n +i > 0}. Then HP 
is a metric space with the hyperbolic distance d : HP x HP — > M. assigning to each (x, y) G HP x M n 
the unique number d(x,y) > such that coshd(x,y) = —(x \ y). The hyperbolic space is uniquely 
geodesic and all its balls are convex. The hyperbolic law of cosines states that in a hyperbolic 
triangle with vertices x,y,z £ HP and 7 the hyperbolic angle between the geodesic segments [x, y] 
and [x, z] we have 

cosh d(y, z) = cosh d(x, y) cosh d(x, z) — sinh d(x, y) sinh d(x, z) cos 7. 

Let i£l and n £ N. The classical model spaces M£ are defined as follows: if k > 0, M^ 1 is 
obtained from the spherical space S n by multiplying the spherical distance with \j\fk\ if k = 0, 
is the n-dimensional Euclidean space E"; and if k < 0, is obtained from the hyperbolic space 
HP by multiplying the hyperbolic distance with l/\/—k. The model spaces inherit their geometrical 
properties from the three Riemaniann manifolds that define them. Thus, if k < 0, MJ} is uniquely 
geodesic, balls are convex and we have a counterpart of the hyperbolic law of cosines. If k > 0, 
there is a unique geodesic segment joining x,y G if and only if d(x,y) < ir/Vk. Moreover, 
closed balls of radius smaller than ■n/{2^/k) are convex and we have a counterpart of the spherical 
law of cosines. We denote the diameter of M" by D K . More precisely, for k > 0, D K = ir/y/Ti and 
for k < 0, D K = 00. 

Now we briefly introduce CAT(A:) spaces. Let (X, d) be a geodesic space. A k-comparison 
triangle for a geodesic triangle A(xi, X21 £3) in (-X^d) is a triangle A(xi, £2, £3) in M| such that 
d M 2(x~i,x~j) = d(xi,Xj) for j G {1,2,3}. For A; fixed, ^-comparison triangles of geodesic triangles 
(having perimeter less than 2-0^) always exist and are unique up to isometry. 

A geodesic triangle A in X is said to satisfy the CAT(k) inequality if, given A a fc-comparison 
triangle for A, for all x,y G A 

d(x,y) < d M 2(x,y), 

k 

where x, y G A are the comparison points of x and y, respectively. 

If A; < 0, a CAT(k) space is a geodesic space for which every geodesic triangle satisfies the 
CAT(/c) inequality. If k > 0, a metric space is called a CAT(/c) space if every two points at distance 
less than can be joined by a geodesic and every geodesic triangle having perimeter less than 
2Dfc satisfies the CAT(/c) inequality. 

Let A be a metric space and C C X. We define the distance of a point z G X to C by 
dist(z,C) = inf d(z,y). The metric projection Pc onto C is the mapping 

Pc{z) = {y G C : d(z, y) = dist(z, C)} for every z £ X. 

The following proposition gathers important properties of the metric projection in the setting 
of CAT(0) spaces. In [13], a counterpart of this result was given in the setting of CAT(A;) spaces 
with k > 0. 
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Proposition 2.1 ([7], Proposition 2.4, p. 176). Let X be a complete CAT(0) space, x G X and 
C C X nonempty closed and convex. Then the following hold: 

1. The metric projection Pc(x) of x onto C is a singleton. 

2. If ye [x,P C (x)], then P c {x) = P c (y). 

3. If x $l C and y £ C with y ^ Pc(x) then Z Pc ^(x,y) > it/2. 

4- The mapping Pq is a nonexpansive retraction from X onto C . Further, the mapping H : 
X x [0, 1] —7- X associating to (x, t) the point at distance td(x, Pc(x)) on the geodesic [x, Pc(x)] 
is a continuous homotopy from the identity map of X to Pq ■ 

We give next a concept of convergence in metric spaces. Let (X, d) be a metric space, (x n ) C X 
a bounded sequence and x £ X. The asymptotic radius of (x n ) is given by 



The sequence (x n ) in X is said to A-converge to x G X if x is the unique asymptotic center of (u n ) 
for every subsequence (u n ) of (x n ). This notion of convergence coincides with the weak convergence 
in Hilbert spaces. For more on this topic and other concepts of weak convergence in metric spaces 
the reader may see [18, 29, 20, 13, 19]. 

3 Main results 

3.1 Reflecting in geodesic spaces 

Let (X, d) be a uniquely geodesic space with the geodesic extension property. Let C be a nonempty 

closed and convex subset of X and suppose the metric projection onto C is well-defined and 

singlevalued. The reflection of a point x G X with respect to C can be any point z in a geodesic 

x ~\~ z 

line containing the geodesic segment [x, Pqx] for which Pqx = — - — . 

If X has no bifurcating geodesies, then geodesies can be extended in a unique way so the 
reflection of a point is uniquely determined. Thus, in such a setting, the reflection mapping Rq 
which assigns to each point its reflection is well-defined and singlevalued. Recall that two geodesies 
bifurcate if they have a common endpoint and coincide on an interval, but one is not an extension 
of the other. Note however that one may also define singlevalued reflection mappings in spaces 
failing this property. 

We study next the reflection mapping in a nonlinear setting. By [16, Example 22.1], we know 
that in the complex Hilbert ball the reflection mapping fails to be nonexpansive which yields that in 
CAT(O) spaces this mapping is not necessarily nonexpansive in general. However, in the particular 
setting of model spaces, the reflection mapping proves to be nonexpansive. This shows that in 
model spaces, as in the case of Hilbert spaces, the metric projection is an averaged mapping with 
constant 1/2. 



r((x n )) = inf < limsup d(x, x n ) : x G X 




and the asymptotic center of (x n ) is the set 
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Proposition 3.1. Let fc£R and n G N. Suppose C is a nonempty closed and convex subset of 
Ml and x,y G M% such that dist(x,C), dist(y,C) < D k /2. Then, 

d(R c x,R c y) < d(x,y). 

Proof. First we consider the case k = — 1. For simplicity, denote c x = Pcx, c y = Pcy, x' = Rqx 
and y' = RcV- Let 7 = Z Cx (x,c y ) and 7' = Z Cy (y,c x ). Notice that, by Proposition 2.1, 7,7' > ir/2. 
Consider the geodesic triangles A(x, c x ,c y ) and A (a;', c x , c y ). By the hyperbolic cosine law we have 
that 

cosh(i(x, Cy) = coshd(x, c x ) cosh d(c x ,c y ) — smhd(x, c x ) sinhd(ca;, c y ) cos 7 

and 

cosh (f(x', Cy) = coshd(x', c x ) cosh d(c x ,c y ) — smhd(x' , c^) sinh d(c x , c y ) cos(7r — 7). 
Since d(x, c x ) = d(x' , c x ) and cos(7r — 7) > 0, we get coshd(x, c y ) > coshd(x', c y ) and thus 

d(x, c y ) > d(x' ', Cy). 

Similarly we get d(y,c x ) > d(y',c x ). 

Consider now the geodesic triangles A(x, c x ,y) and A(x', c x , and denote f3 = ^ Cx (x, y). Applying 
again the hyperbolic cosine law we obtain that 

coshd(x,y) = cosh d(x, c x ) cosh d(c x , y) — sinh d(x, c x ) sinh d(c x , y) cos /3 

and 

cosh <i(x', y) = cosh d(x' , c^) coshd(c x , y) — sinh <i(x', c x ) sinh d(c x , y) cos(7r — [3). 
Adding these two equalities we get that 

cosh d(x, y) + cosh d(x',y) = 2 cosh d(c x , x) cosh d(c x , y). (3-1) 

In a similar way we have that 

cosh d(x' , y') + cosh d(x, y') = 2 cosh d(c x , x) cosh d(c x ,y'), (3-2) 

cosh <i(x, y) + coshd(x, y') = 2 cosh d(c y , x) cosh (/(c^, y), (3-3) 
cosh d(x', y') + cosh d(x' , y) = 2 cosh d(c y ,y) cosh d(c y ,x'). (3.4) 
Suppose now that d(x',y') > d(x,y). From (3.1) and (3.2), 

coshd(x', y') + coshd(x,y') = 2coshd(c x ,x) coshd(c x ,y') 

< 2 cosh d(c x , x) cosh d{c x , y) 

< cosh ci(x', y') + cosh d(x', y), 

which implies that d(x,y') < d(x',y). By using (3.3) and (3.4) in a similar way, we get that 
d(x,y') > d(x',y), which is a contradiction and thus the result follows. The proof in any other 
model space MJ} follows similar patterns by using the corresponding law of cosines of each space. □ 

The above result also holds in the equivalent infinite dimensional spherical and hyperbolic spaces 
defined using points in £2 (see, for instance, M.°° in [14] and [13, Example 3.4]). 

In the sequel, we analyze the reflection mapping in some particular CAT(/c) spaces. More 
precisely, we consider the gluing of two model spaces M£ and MfJ (for more details on gluings see 
[7, Chapter 1.5]). Note that two spaces of constant curvature can only be glued through geodesic 
lines, geodesic segments or singletons. First we consider the gluing of two model spaces by a 
singleton. Although the reflection mapping is not uniquely defined in this type of gluing spaces, 
there exists a natural way of defining it such that it is nonexpansive. 
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Proposition 3.2. Let k,k' G M., k ^ k' and m,n G N. Consider the gluing space (X,d) of MJ} 
and MJ}} through a point 9. Suppose B is a nonempty closed and convex subset of X and x G X 
such that dist(x,B) < min{Dfc, Dk'}/2. We define the reflection of a point with respect to B in X 
in the following way: 

(1) If 9 G B , Rbx = R BnM px, where M? G {M£,MJ}}} is the model space that contains x and 
Rbdm v i s th e reflection in the model space Mj with respect to B n Mj . 

(2) If 9 ^ B, observe that B is strictly contained in just one of the two model spaces that define 
the gluing X. Suppose B C MJ}. Then, 

(2.1) if x G MJ}, Rbx is the point obtained by reflecting in MJ} with respect to B; 

(2.2) if x G M™, the reflection is uniquely defined since d{9,Psd) > 0. Specifically, Rbx is the 
point in the geodesic line in MJ} containing [9,P B 9] such that d{Ps9, Rbx) = d(x,PB9) 
and d(R B x, 9) > d{Rpx, Pb&)- 

Then Rb is nonexpansive. 

Proof. (1) Suppose 9 G B. We prove that Rb defined as in (1) is nonexpansive. Let x,y G X. The 
cases x,y G MJ} or x,y G M^} are immediate since Rb reduces to the reflection in model spaces. 
Suppose now x G MJ} and y G M™. Then, Rbx = R B r\M%x and Rsy = RbciM^U- Consider the 
geodesic triangles Ai = A(Rbx,Pbx,9) and A2 = A(x, Pbx,9). Notice that Ai,A2 C MJ}. For 
simplicity, we consider k = — 1 (the proof for k 7^ — 1 follows similar patterns). Since Zp BX (x,9) > 
7r/2, similarly as in the proof of Proposition 3.1, d(9,x) > d(9, Rbx). Likewise, we can see that 
d(9,y) > d{9,Rsy)- Therefore, 

d(R B x, R B y) = d{R B x, 9) + d{0, R B y) < d(x, 9) + d{y, 9) = d(x, y), 

and the conclusion follows. 

(2) Next we suppose that 9 ^ B and B C M£. We prove that Rb defined as in (2) is nonexpansive. 
Let x,y G X. As we mentioned before, the case x,y G MJ} is immediate. Let x, y G MJ}}. Since 
dist(0, B) > 0, we have that Pbx = Psy = Pb& and so 

d(R B x, R B y) = \d(x, P B x) - d(y, P B y)\ = \d(x, 9) - d(y, 0)\ < d(x, y). 

Finally, suppose x G MJ}, y G MJ}}. As in the previous case, Psy = PbQ- Consider z = Rp B oRsy G 
MJ}, where Rp B g is the reflection mapping with respect to the point Pb9 inside MJ}. By definition 
we have Pb9 G [RbUi Since in a model space geodesies do not bifurcate, the geodesic line 
containing [R B y, P B 9] coincides with the geodesic line containing [R B y, 9]. Let c be this geodesic 
line. Then, by definition of the reflection, z G c. Moreover, d{Rsy, Pb&) = d{z,Ps9) > d{Ps9,9) 
and d{9,Rsy) > d{Rsy, Pb@)- Thus, 9 G [Pb9,z]. As a consequence, 

d(P B 9,y) = d{P B 9,R B y) = d(P B 9,z) = d(P B 9 , 9) + d{9 , z) . 

On the other hand, P B z = Pb9 and therefore Rbz = RbV- Finally, since d(x, z) < d(x, 9)+d(z, 9) = 
d(x, y) — d(y, 9) + d(z, 9) = d(x, y) and the reflection in MJ} is nonexpansive, we get 

d{Rsx, Rsy) = d{Rsx, Rbz) < d{x, z) < d(x, y), 

and the result follows. □ 
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Next we see that when gluing two model spaces through a geodesic segment or a geodesic line, 
we cannot define the reflection in the gluing in such a way that it is nonexpansive with respect to 
every nonempty closed and convex set. To illustrate this fact we consider the particular gluing of 
H 2 and R 2 . 

Example 3.3. Let c: R — > M 2 be the geodesic line in the plane defined by c(t) = (t, 0) and c' : K — > 
H 2 the geodesic line in the hyperbolic plane defined by c'(t) = (0, sinht, cosht). We consider the 
gluing of these two model spaces by the correspondence given by the isometry j : c(R) — > c'(K) 
defined as j(c(t)) = c'(t) for every t £ R. 

Let x = (— arccosh\/2, h) and y = (arccosh\/2, h), where h > 0. Take p = (\/ z 2 — 1,0, z), with 
z > 1 and B = {p}. 

We see that there exist values of h and z such that d(R^x, Rbu) > d(x,y) = 2arccosh\/2. Notice 
that 

d(y,p) = mi\d(y,c(t)) + d(c'(t), p)\ = inf < yfarccoshv^ — t) 2 + h 2 + arccoshfz cosh t) >. 

iGR teK [ v J 

Let f(t) = d(y,c(t)) +d(c'(t),p). Then, 

z sinh t arccosh\/2 — t 



V z 2 cosh 2 1 - 1 J (arccoshv 7 ^ - t) 2 + h 2 



f'(t) 



Since lim / (t) = 2, lim f'(t) = —2 and f is increasing in t, then there exists only one value 

t— >oo t— >— oo 

to £ K that gives the infimum. Moreover, since f'(0) < and /'(arccosh-v/2) > 0, then to £ 
(0, arccosh\/2) . 

By using the symmetry of the space, we see that if to is the value where the previous infimum is 
attained, then —to is the value that gives d(x,p). Moreover, d(x,p) = d(y,p). 

Let d = d{Rsx, Rbv) o,nd 7 be the hyperbolic angle between the segments \p, c(to)] and [p, c(— to)]. 
By the hyperbolic cosine law, we have 

cosh 2to = (cosh d(p,c (to))) 2 — (sinh d(p, c(to))) 2 cos 7 

and 

cosh a! = (cosh d(x,p)) 2 — (sinh d(x,p)) 2 cos 7. 

If we consider for instance h = 1/100 and 1 < z = (\^2 + 15)/16 < then to — 0.8392 and 
cos 7 ~ —0.7991. Consequently, coshd ~ 3.7363 > cosh d(x,y) = 3. 

The underlying idea in this construction is the fact that geodesies in M 2 diverge faster than in M. 2 . 

Since the previous construction can be easily adapted to the case of gluings through segments, 
we also obtain the same conclusion for this type of gluings. 

3.2 Convergence results 

We begin this section by defining the averaged alternating reflection (AAR) method in the setting 
of geodesic spaces. Let (X, d) be a uniquely geodesic space. Given A and B two nonempty closed 
and convex subsets of X, suppose that the reflection mappings Ra and Rb are well-defined and 

singlevalued. Consider the mapping T: X — > X defined by T = — ^ — — . The AAR method 
generates the following sequence for a starting point xq £ X: x n = T n xo for every n > 1. 
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We prove next the convergence of the AAR method in the setting of model spaces with k < 
which is an analogue of a weak convergence result in Hilbert spaces [4, Fact 5.9]. Note that we 
obtain strong convergence since the model spaces are proper metric spaces (see also [19, Proposition 
4.4]; recall that a space is proper if every closed ball is compact). 

Theorem 3.4. Let k < and n G N. Suppose A and B are two nonempty closed and convex 
subsets of MJ} with Af\ B ^ 0. Let xq G M£ and (x n ) be the sequence starting at xq generated by 
the AAR method. Then 

(1) (x n ) converges to some fixed point x of the mapping T and Pbx G A n B. 

(2) The "shadow" sequence (Pex n ) is convergent and its limit belongs to A(~)B. 

Proof. (1) By Proposition 3.1, RaRb is nonexpansive which yields that T is also nonexpansive 
because M£ with k < is Busemann convex. In addition, Fix(T) ^ since A n B ^ 0. Thus, 
all orbits of T are bounded. Using [15, Proposition 2] it follows that T is asymptotically regular. 
Applying [1, Proposition 6.3], we get that (x n ) A— converges to a fixed point of T, which implies 
that (x n ) converges to a fixed point of T. Let x be the limit of (x n ). Since x = Tx it follows that 

x = RaRbx. Moreover, because Pbx = , 

RaRbx + Rbx 
P B x = = PaRbx g A. 

Hence, P B x G A n B. 

(2) is immediate since the metric projection Pg is continuous in CAT(O) spaces. □ 

We remark that a similar result holds in the gluing spaces described in Proposition 3.2 (when 
k, k' < 0). In H°° one obtains an analogous A-convergence result. 

In model spaces MJ} with k > the mapping T used in defining the AAR method may fail to 
be nonexpansive as the following example shows. 

Example 3.5. Consider the points a = (1,0,0), b = (v2/2, v2/2, 0) and p = (0,0,1) in the 
spherical space S 2 . Let c G [a,p] and d G [b,p] such that d(a,c) = d(b,d) = tt/8. Suppose also 
A = [a, b] and B = [c, d] . To see that the mapping T is not nonexpansive it is enough to note that 
Tc = P A c, Td = P A d, but d(Tc, Td) = d{P A c, P A d) = d(a, b) = vr/4 > d(c, d). 
Taking B to be the whole positive octant, we obtain a similar example for the case Af] B ^ 0. 

However, for model spaces MJ} with k > we can prove the result below. Note that the 
existence of a convergent subsequence of (x n ) already guarantees that A n B ^ 0. 

Theorem 3.6. Let k > and n G N. Suppose A and B are two nonempty closed and convex 
subsets of M£. Let C C MJ} be nonempty convex with diam(C) < Dk/2 such that A,BOC and 
RaRb{C) C C . Let xq G C and (x n ) be the sequence starting at xq generated by the AAR method. 
Then, 

(1) Any convergent subsequence (x nfe ) of (x n ) converges to a fixed point x of the mapping T and 
P B x G An B. 

(2) The "shadow" sequence (PsXn) is bounded and each of its cluster points belongs to All B. 
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Proof. (1) Since diam(C) < D/./2, we know that the metric is convex on C. Reasoning as in the 
previous theorem, we get that d(x n ,x n+ i) — > as n — > oo. Let (x nk ) be a convergent subsequence 
of (x n ) and x be its limit. Since RaRb is continuous we have RaRbx^ — > RaRbx. Besides, 
(x nfe +i) also tends to x, which implies that x = Tx. The fact that Pbx € An B follows as in the 
proof of the previous theorem. 

(2) Let (PBX nk ) be a convergent subsequence of (PsXri) an d y be its limit. We may suppose that 
(x n , k ) is convergent (otherwise consider a convergent subsequence of it) and denote x = lim x nk . 

k— >oo 

From the continuity of the metric projection Pg, we get that y = Pbx and we are done. □ 

A more general way of defining the sequence (x n ) would be to use the Krasnoselski-Mann [26, 23] 
iteration starting at xq £ X: 



-1 = (1 - X n )x n + X n RARBX n , 



(3.5) 



where (A n ) C [0,1]. 

Results similar to Theorems 3.4, 3.6 hold for the sequence (x n ) generated by (3.5) when assuming 
for instance that (A n ) is divergent in sum and bounded away by 1. In this case, consider F = 
Fyx.{RaRb) 7^. AC\B and note that (d(x n ,p)) is decreasing for each p £ F. Apply [15, Proposition 
2] as above to get that d(x n , RaRb x u) — > (if k < 0, use asymptotic center techniques similar to 
those considered in [1, Proposition 6.3] to obtain the convergence of (x n )). Notice that now one 
obtains convergence to a fixed point of the mapping RaRb (see also [12]). 

We finish this section by giving a rate of asymptotic regularity for the sequence (x n ) generated 
by (3.5). Using proof mining methods, Kohlenbach [21] and later Kohlenbach and Leu§tean [22] 
computed exponential (in 1/e) rates of asymptotic regularity for the Krasnoselski-Mann iteration 
in normed and hyperbolic spaces, respectively. The next result which follows from [22, Corollary 
3.18] gives an explicit bound on the rate of asymptotic regularity for the sequence (x n ). 

Theorem 3.7. Let k £ R and n £ N. Suppose A and B are two nonempty closed and convex 
subsets of MJ} with A n B ^ 0. Let C C MJ} be nonempty convex and of diameter bounded 
above by b G (0,D k /2) such that A,B<ZC and RaRb{C) C C. Assume K £ N, K > 2 and 

Then, 



(X n ) C 



1 

K' 1 



K 



Vx G X,Ve > 0,Vn > $(e,b,K), d(x n ,x n+ i) < e, 



where 



with 



<S>(e,b,K) := KM 



M > 



(K — !)(! + 26) 
Ke 



Thus, the above rate only depends on e, on an upper bound b on the diameter of the set C and on 
(A n ) via K. Rates of asymptotic regularity for the Krasnoselski-Mann iteration were further studied 
and improved by Leu§tean [24] in the setting of uniformly convex hyperbolic spaces that admit a 
modulus of uniform convexity which decreases with respect to the radius. As a consequence, one 
gets a quadratic in 1/e rate of asymptotic regularity in CAT(O) spaces for constant X n = A £ (0, 1) 
[24, Corollary 19]. In particular, these results can be applied in the setting of model spaces MJ} 
with k < 0. 
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4 Appendix: A related property of model spaces 

We have seen above that the reflection mapping with respect to nonempty closed and convex subsets 
is nonexpansive in the context of model spaces. The property we prove below implies in particular 
that the reflection mapping is nonexpansive in all model spaces M£ with k < 0. 

Proposition 4.1. Let k G R, k ^ and n G N. Consider C C MJ? with diam(C) < and the 
points x, y, x' , y',a, b G C for which 

a= -x + -x and b = -y + -y . 

Denote byy andy' the angles between the geodesic segments [y,b] and [b,a], [x,a] and [a,b\. Suppose 
that 

d(a, b) < d ((1 - X)x + Xa, (1 - X)y + Xb) , for each X G [0,1]. (4.6) 
If one of the following holds 

(i) d(x, a) = 0, 

(ii) d(x,a) = d(y,b), 

(hi) d(x, a) < d(y, b) and (j < n/2 or i > tt/2 for k < 0), (y > tt/2 or j' < tt/2 for k > 0), 
then d(x',y') < d(x,y). 

Proof. We prove the result for H n , the proof for the other model spaces can be obtained by using 
the corresponding law of cosines. 
For A G (0, 1), consider 

x\ = (1 — X)x + Xa and y\ = (1 — X)y + Xb. 

Denote d(x, a) = A and d(y, b) = B. 

(i) : By the hyperbolic cosine law, 

coshd(a,6) < cosh<i(a,yA) = coshd(a, b) cosh d(y\, b) — sinhd(a, b) sinh d(y\, b) cos 7. 

Thus, 

•1^ K\ CQ8h((l-A)^)-l 

sinh ata, b) cos 7 < —-. — cosh a(a, 0). 

sinh ((1 — X)B) 

Letting A — > 1 we obtain that cos 7 < 0. Using again the hyperbolic cosine law it follows that 
coshd(a,y') < coshd(a,y), hence, d(a,y') < d(a,y). 

(ii) or (iii): Apply again the hyperbolic cosine law in the triangles A(a,y,b) and A(a,y',b) to 
obtain that 

cosh d(a, y) + cosh d(a, y') = 2 cosh d(a, b) cosh B. (4.7) 
By the hyperbolic cosine law applied in the triangles A(x\,y,y\) and A(x\,y' ,y\), 

sinh ((2 — X)B) coshd(x\,y) + sinh (XB) coshd(x\,y') = sinh (2B) coshd(x\,y\). (4-8) 

Using again the hyperbolic cosine law in the triangles A(a,x,y) and A(x\,x,y), 

cosh d(a, y) — cosh d(x, y) cosh A cosh d(x\,y) — cosh d(x, y) cosh (XA) 
sinh A sinh (A A) 
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from where 

sinh (XA) cosh d(a, y) + sinh ((1 — A) A) coshd(x, y) = sinh A cosh d(x\,y). 
Applying a same reasoning in the triangles A(a, x',y') and A(x\,x' ,y') we obtain 

sinh ((2 — X)A) cosh d(a, y') = sinh A cosh d(x \, y') + sinh ((1 — X)A) cosh d(x', y'). 
This implies 

sinh (XA) sinh ((2 — X)B) (cosh d(a, y) + cosh d(a, y')) 

+ cosh d(a, y') (sinh ((2 - X)A) sinh (AS) - sinh (XA) sinh ((2 - X)B)) 
+ sinh ((1 - X)A) sinh ((2 - A)5) cosh d(x, y) 
= sinh A (sinh ((2 — A)B) cosh.d(x\,y) + sinh (AS) cosh<i(xA, y')) 
+ sinh ((1 - X)A) sinh (AJ3) coshd(x', ?/')• 

Using (4.7) and (4.8) and the fact that cosh d(x\, y\) > coshd(a,b) we obtain 

2 cosh B cosh d(a, b) (sinh (XA) sinh ((2 - X)B) - sinh A sinh B) 

+ cosh d(a, y') (sinh ((2 - A)A) sinh (XB) - sinh (A4) sinh ((2 - X)B)) 
+ sinh ((1 - X)A) sinh ((2 - X)B) cosh y) 

> sinh((l - A) -A) sinh (XB) cosh j/). 

Dividing by sinh((l — A)A) and letting A — > 1, 

2 (cosh £> cosh d(a, b) — cos d(a, y')) ^— sinh yl cosh B — cosh A sinh B^j 

> sinhi? (coshd(x ! , y') — coshd(x,y)) . 

Therefore, 

B 



2 sinh d(a, b) cos(7r — 7) ^— sinh^lcoshi? — coshAsinhSj > cosh d(x', y') — cosh d(x, y). 
In a similar way we obtain that 

2 sinh(i(a, b) cos(7r — 7') ^— sinh B cosh j4 — cosh B sinh ^4^ > cosh d(x' , y') — cosh d(x, y). 

Since x 1— > is non-increasing, the conclusion follows. 



□ 



Note that when k = the above property is immediate. In fact, in any Hilbert space, it is easy 
to see that condition (4.6) already implies that d(x',y') < d(x,y) (see also the so-called property 
(S) with 6 = 1 studied in [8] which is equivalent to this condition in the setting of normed spaces). 
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